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We study the zeros of the partition function in the complex j8 plane (Fisher's zeros) in SU(2) 
and SU(3) gluodynamics. We discuss their effects on the asymptotic behavior of the perturbative 
series for the average plaquette. We present new methods to infer the existence of these zeros in 
region of the complex j3 plane where MC reweighting is not reliable. These methods are based on 
the assumption that the plaquette distribution can be approximated by a 0^ type distribution. We 
give new estimates of the locations for a 4^ lattice. For SU{2), we found zeros at j3 = 2.18(1) ± 
/0.18(2) (which differs from previous estimates), and at j3 = 2.18(1) ± /0.22(2). For SU{3), we 
confirm J3 =5.54(2)±i0.10(2)and found additional zeros at J3 =5.54(2)±/0.16(2). Some of the 
technical material can be found in recent preprints, in the following we emphasize the motivations 
(why it is important to know the locations of the zeros) and the challenges (why it is difficult to 
locate the zeros when the volume increases). 
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1. Introduction 

More than 200 years ago, L. Euler realized that the asymptotic behavior of the expansions 
coefficients of the Bessel functions about the origin could be obtained from the location of the 
complex zeros. For a single plaquette with SU (2) gauge group (as defined, for instance, in Ref. 
[^, as usual j3 =4/g^), we obtain that in this simple model the average plaquette P can be written 
as 

where the a„ are the locations of the zeros of J\ {£) (on the positive real axis). Using the geometrical 
series one sees that the radius of convergence of the expansion in j3 is 1 /tti . This result can be 
rephrased in the following way: the radius of convergence of the strong coupling expansion of P 
is the distance from the origin to the closest zero of the partition function (which in this simple 
example is located, as well as all the other zeros, on the imaginary axis). 

The situation is more interesting at large j8 and fundamentally different from the case of scalar 
models. As well known, the zeros of J\ (z) become approximately equally spaced for large real 
argument. If we now consider P in the = 4/j8 complex plane, the zeros accumulate at zero 
along the imaginary axis which plays the role of a Stokes line. Unlike the scalar case the partition 
function is well defined at real negative g^, but as we cross the imaginary axis, the values of P 
changes discontinuously. This can be seen from the sum rule P{^) +P(— j3) = 2, which remains 
true on arbitrary lattices with even number of sites in every directions [^]. 

We expect that the features of the strong and weak coupling expansions observed for the one 
plaquette model to remain true on generic lattices. However, other zeros of the partition function 
related to phase transitions are expected. For the 2-dimensional Ising model, it has been observed 
by M. Fisher [Q] that the zeros of the partition function pinch the real axis of the complex inverse 
temperature at jS^ . For the 3-dimensional Ising model, numerical calculations consistent with this 
scenario have been obtained [Q] on a 4^ lattice. These zeros in the complex j8 plane are called 
Fisher's zeros and should be distinguished from the Lee- Yang zeros in the complex magnetic field 
plane. 

For a zero temperature lattice gauge theory with a Wilson action, we expect no phase transition 
on the real j8 axis. However, such theories can be seen as "close" to other theories (e. g. at non-zero 
temperature or with a positive adjoint coupling) that have a phase transition. Consequently, it is 
plausible that zero temperature lattice gauge theory with a Wilson action have Fisher zeros close 
to the real axis but that these zeros do not pinch the real axis in the infinite volume limit. This is 
a plausible explanation for the unexpected behavior of the weak coupling expansion of P 
for SU (3) [||, ^]. Standard methods of series analysis suggest ^] a singularity on the real axis, 
namely P oc (1/5.74 — l/jS)^"^. This would imply a peak in the second derivative of P with a 
height increasing with the volume, which is not seen at zero temperature [^. The vicinity of the 
critical point in the fundamental-adjoint plane, suggests the approximate mean field behavior [^]: 

-5/'/5i3ocin((l/j3„-l/j3)2 + r2) , (1.2) 

Fits of the series with such parametric form yield the approximate values j3,„ ~ 5.78 and Y ~ 0.006 
(i.e Im j8 ~ 0.2). 
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Values of T which are too large (too small) would produce modulations of the coefficients 
(peaks in the derivatives of P) which are not observed. A detailed analysis [§] yield the bounds (for 
SU (3)) 0.001 < r < 0.01. This suggests zeroes of the partition function in the complex j3 plane 
with 

0.03 ~ O.OOljS,^, < /mjS < 0.01j3^ ~ 0.33 (1.3) 

In these proceedings, we report recent efforts to locate such zeros in pure gauge theories with 
gauge group SU (2) and SU (3). Some of the technical details can be found in a recent preprint [pl|]. 



Another preprint is in preparation [12] and should be available soon. In order to avoid repetitions, 



we wiU emphasize the motivations and aspects not covered in these preprints. 

2. The non-perturbative part of the plaquette and the gluon condensate 

It is quite common that the difference between a physical quantity and its perturbative ex- 
pansion is of the form exp(— /T/g^). One of the best known example is the quantum mechanical 
double-well where the perturbative series is not able to take into account the tunneling effect and 
instantons are needed. For the average plaquette, the issue is obscured by the hypothetical zero 



close to the real axis and the factorial growth necessary to get an envelope [10] in the accuracy 
versus coupling at successive order is not reached at the order where the perturbative expansion is 
available. Larger order extrapolation are necessary. Two models have been considered [p^]. 
The first is based on Eq. ( |l.2[ ) which implies a dilogarithmic series for P: 



£ au^-'' ~ C(Li2(j8- V(i8™ ' + /r)) + h.c . (2.1) 



k=0 



The low order coefficients depend very little on Y as long as the upper bound given above is sat- 
isfied. It is remarkable that the method of stochastic perturbation theory allows us to calculate the 
series up to order 10 [||] and 16 [^, however larger series are needed to resolve F! As we will ex- 
plain below, the situation might be better for SU (2). The two other parameters can be determined 



using for instance the values of ag and aio. This gives very good predictions [ ]10| ] of the values of 
a^,a-i,.. .. This nice regularity is still begging for a diagrammatic explanation. 

Despite its predictive success, the dilogarithmic series has a finite radius of convergence and 
the coefficients do not have the expected factorial growth that is observed in the one plaquette 
model. The second extrapolation was based on an IR renormalon model [|lj, 14,15] 

Yatt'^K Fdlc-'^' (l-?33/167r2)-'-204/i2i ^2.2) 

~^=^{\+di/^ + ...) (2.3) 
These two extrapolations seem consistent with the behavior 

P{^)-Ppertm~C{a/rof (2.4) 

with a(j8) defined with the force scale [[T^, |l^] with ro = 0.5 fm, and Ppen appropriately truncated. 
For large j8 this has the desired exponential form. Attempts have been made in the past [|T8], ^ ^ 
to relate C to the so-called gluon condensate JT4l]. Remembering the a/n factor in the definition. 
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Figure 1: Distribution of the 500,000 values of S in an histogram with 100 bins for a SU{2) pure gauge 
theory on a 4'* lattice at j3 = 2.20. The solid line is a Gaussian fit. 

the value that could in principle be compared with the commonly used value of 0.012 GeV^ is 
(36/71^ )CrQ'* for Nc = 3. C is sensitive to resummation. C ~ 0.6 with the bare series [10] and 0.4 



with the tadpole improved series [Q]. This gives values 3-5 times larger than the value quoted above. 
Besides the question of scheme dependence, the gluon condensate is not an order parameter and it 
seems difficult to compare the lattice results with quantities defined in the context of sum rules. On 
the other hand, it is important to figure out how well the scaling with the lattice spacing given in 



Eq. (2.4) is obeyed and if it can be explained semi-classically as it can be done for the double- well. 
In this context, it is crucial to understand the complex singularities of P, which complicate the 
analysis of the scaling. 

3. The zeros of the partition function 



The zeros of the partition function can be located with the reweighting method [19, 20 1 



Z(j8o + Aj8)=Z(i3o)<exp(-Aj85)>ft . (3.1) 

It is convenient to subtract < S > from S in the exponential because it removes fast oscillations 
without changing the complex zeroes. Z(j8) is the Laplace transform of density of states n{S): 

poo 

Z(j8)= / dSn{S)&x^{-^S) . (3.2) 

JO 

If we could estimate Z(j8) along an axis in the imaginary direction, we could calculate the inverse 
Laplace transform and obtain n{S). As shown in Fig. |l|, the distribution of values of S for SU (2) 
can be fitted very well with a Gaussian with a| =< > — <S >^. This suggests [20] a criterion 
to determine a region of confidence for MC zeros. For a Gaussian distribution, the fluctuation in 
exp(— Aj3(5— < S >)) is smaller than the average for |Aj8p < \n{Nconf)/<ys- "^^^^ defines aradius 
of confidence y^ln{Nconf.)/ <^s in the complex j3 plane. As a| oc V, the radius of confidence shrinks 
like V^^/^ which is hard to beat with y^lnNconf.- 

When the zeros pinch the real axis as the volume increases, which is expected near a phase 
transition, there is some hope that it is possible to follow them with MC methods. This explains 
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Figure 2: n = [Ni-NPi)/ ^/Wi for SU{2) at j3 = 2.20 on 4^ and 6"^ lattices. 



some positive results pO| ] on 2 x and 4 x L^* lattices. On the other hand, on lattices, we 
expect the zeros to stabilize away from the real axis and MC methods can only give lower bounds 
that shrink when L increases. For instance for SU (3) on a 8^ lattice [21], we were only able to 



verify that /mj3 > 0.03. In addition, for SU (2), there seem to be no zeros in the Gaussian region of 
confidence even on a 4^ lattice. For these reasons we have developed new methods to find zeros of 
the partition that lay outside of the region of confidence of MC calculations. 

4. New methods to locate the zeros |l^] 

Gaussian distributions (of S) have no complex zeros. The Gaussian circle of confidence in the 
complex j8 plane defined by the condition Of < \ f\ sjNconf.- If this criterion is applied directly to 
a non-Gaussian distribution having complex zeros, it will automatically exclude the regions that 
contain the zeros. When looking for complex zeros, we look for the intersection of the zero level 
curves for the real and imaginary parts. We are interested in knowing how much these level curves 
can move due to statistical fluctuations. We proposed JTl] ] to consider the alternative region of 
confidence defined by a condition that controls the error on the level curves: 



a/ < ^/y^W l/'l • (4-1) 
In order to be useful d should be a fraction of the typical distance between zero level curves of the 



real and imaginary part. This methods has allowed [11] to reject dubious zeros on the edge of the 
Gaussian circle of confidence for SU (2) on a 4^ lattice. It remains applicable when the deviation 
from a Gaussian distribution is significant and true zeros appear. 

The zeros come from the deviations from the Gaussian behavior. As shown on Fig. discrep- 
ancies in unit of the expected fluctuations are coherent for L = 4 but as the the volume increases, the 
signal gets lost in the noise (for that particular value of j8). The nice regularities of the difference 
with the Gaussian approximation (for small lattices) suggest to fit the distribution with 

P{S) oc exp(-Ai5 - A25^ - hS^ - hS'^) (4.2) 

The unknown parameters were determined from the first four moments using Newton's methods 
and also by minimization. Very good agreement between the two methods was found on 4^ 
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Figure 3: Zeros of the real (crosses) and imaginary (circles) using MC on a 4^ lattice, for SU{2) at j3 = 
2.18 and SU{3) at j3 =5.54. The smaller dots are the values for the real (green) and imaginary (blue) parts 
obtained from the 4 parameter model. The MC exclusion region boundary for li = 0.15 is represented by 
boxes (red). 



lattices. The zeros can then be calculated for the parametric form (4.2) using accurate numerical 
integration. The results are shown in Fig. |3| on a 4^ lattice. Comparing results at different j3 on a 
4^ lattice, we obtained the following locations of the complex zeros: 



j3 = 2.18(1) ±/0.18(2) forSU{2) (differs from [19] 2.23±/0.155 obtained withMC outside 



regions of confidence) and another zero at j8 = 2.18(1) it /0.22(2). 



• j8 = 5.54(2) ±/0. 10(2) for SU{3) (agrees with ^) and another zero at j3 = 5.54(2) ± 
/0.16(2) . 

Note that the ratio of the imaginary and real parts of the closest zero is almost 5 times larger in 
SU (2). This indicates that modulations in the perturbative coefficients of P should be easier to see 
than in SU{3). 



5. Conclusions 

We have build a ladder of methods that can be applied for increasing values of the imaginary 
part. We found new ways to distinguish fake and true MC zeros that work well with non-Gaussian 
examples. Fitting methods based on cubic and quartic perturbations give consistent results at dif- 
ferent P for larger values of the imaginary part on a 4^ lattice. Results on larger lattices will be 
available soon We are in the process of checking the selfconsistency of the parametriza- 

tion at different j8 and are attempting to extract the density of states. Effect of an adjoint term, 
finite-temperature and decimation are also under study. 
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